Ordering of Elements for the Volume & Neighbors Algorithm Constructing Elimination Trees for 2D and 3D h-adaptive FEM

AbstractIn this paper we analyze the optimality of the volume and neighbors algorithm constructing elimination trees for three dimensional h-adaptive finite element method codes. The algorithm is a greedy algorithm that constructs the elimination trees based on the bottom up analysis of the computational mesh. We compare the results of the volume and neighbors greedy algorithm with the global dynamic programming optimization performed on a class of elimination trees. The comparison is based on the Directed Acyclic Graph (DAG) constructed for model grids. We construct DAGs for two model grids: a two dimensional grid refined towards point singularitiy and a two dimensional grid refined towards edge singularity. We show that the quasi-optimal trees created by the volume and neighbors algorithm for the model grids are also captured by the dynamic programming procedure. It means that created elimination trees are optimal in the considered class of elimination trees. We show that different element orderings at the input of the volume and neighbors algorithm result in different computational costs of the multi-frontal solver algorithm executed over the resulting elimination trees. Finally we present the ordering of elements that results in optimal (in the considered class) elimination trees. The theoretical results are verified with numerical experiments performed on a three dimensional grids with point, edge and face singularities

Application of a Hierarchical Chromosome Based Genetic Algorithm to the Problem of Finding Optimal Initial Meshes for the Self-Adaptive hp-FEM

The paper presents an algorithm for finding the optimal initial mesh for the self-adaptive hp Finite Element Method (hp-FEM) calculations. We propose the application of the hierarchical chromosome based genetic algorithm for optimal selection of the initial mesh. The selection of the optimal initial mesh will optimize the convergence rate of the numerical error of the solution over the sequence of meshes generated by the self-adaptive hp-FEM. This is especially true in the case when material data are selected as a result of some stochastic algorithm and it is not possible to design optimal initial mesh by hand. The algorithm has been tested on the non-stationary mass transport problem modeling phase transition phenomenon

Hypergraph Grammars in hp-adaptive Finite Element Method

AbstractThe paper presents the hypergraph grammar for modelling the hp-adaptive finite element method algorithm with rectangular elements. The finite element mesh is represented by a hypergraph. All mesh transformations are modelled by means of hypergraph grammar rules. These rules allow to generate the initial mesh, to assign values of polynomial order to the element nodes, to generate the matrix for each element, to solve the problem and to perform the hp-adaptation

Hypergrammar-based parallel multi-frontal solver for grids with point singularities

This paper describes the application of hypergraph grammars to drive linear computationalcost solver for grids with point singularities. Such graph grammar productions are the rstmathematical formalism used to describe solver algorithm and each of them indicates thesmallest atomic task that can be executed in parallel, which is very useful in case of parallelexecution. In particular the partial order of execution of graph grammar productions can befound, and the sets of independent graph grammar productions can be localized. They canbe scheduled set by set into shared memory parallel machine. The graph grammar basedsolver has been implemented with NIVIDIA CUDA for GPU. Graph grammar productionsare accompanied by numerical results for 2D case. We show that our graph grammar basedsolver with GPU accelerator is order of magnitude faster than state of the art MUMPSsolver

Petri Nets Modeling of Dead-End Refinement Problems in a 3D Anisotropic hp-Adaptive Finite Element Method

We consider two graph grammar based Petri nets models for anisotropic refinements of three dimensional hexahedral grids. The first one detects possible dead-end problems during the graph grammar based anisotropic refinements of the mesh. The second one employs an enhanced graph grammar model that is actually dead-end free. We apply the resulting algorithm to the simulation of resistivity logging measurements for estimating the location of underground oil and/or gas formations. The graph grammar based Petri net models allow to fix the self-adaptive mesh refinement algorithm and finish the adaptive computations with the required accuracy needed by the numerical solution

Applications of a hyper-graph grammar system in adaptive finite-element computations

This paper describes application of a hyper-graph grammar system for modeling a three-dimensional adaptive finite element method. The hyper-graph grammar approach allows obtaining a linear computational cost of adaptive mesh transformations and computations performed over refined meshes. The computations are done by a hyper-graph grammar driven algorithm applicable to three-dimensional problems. For the case of typical refinements performed towards a point or an edge, the algorithm yields linear computational cost with respect to the mesh nodes for its sequential execution and logarithmic cost for its parallel execution. Such hyper-graph grammar productions are the mathematical formalism used to describe the computational algorithm implementing the finite element method. Each production indicates the smallest atomic task that can be executed concurrently. The mesh transformations and computations by using the hyper-graph grammar-based approach have been tested in the GALOIS environment. We conclude the paper with some numerical results performed on a shared-memory Linux cluster node, for the case of three-dimensional computational meshes refined towards a point, an edge and a face

Non-carious tooth loss in terms of erosion – a literature review

The loss of tooth structure caused by the process of chemical dissolution and mechanical abrasion is becoming a great problem in modern dentistry. Patient awareness as well as diagnostic alertness of a dentist are crucial for proper prevention and treatment. Late diagnosis and treatment of extensive loss of tooth structure caused by erosion, abrasion or abfraction poses many difficulties for clinicians. Chronically progressive loss of tooth structure can cause loss of vertical dimension of the bite, secondary orthodontic complications and also aesthetic and functional problems. The treatment is often associated with long-term diagnosis and interdisciplinary, expensive rehabilitation of the entire stomatognathic system. The aim of this study was to review current publications on the prevalence of tooth wear in developmental age, with focus on etiology, risk factors, prevention and treatment methods

Heuristic algorithm to predict the location of C^{0} separators for efficient isogeometric analysis simulations with direct solvers

We focus on two and three-dimensional isogeometric finite element method computations with tensor product Ck B-spline basis functions. We consider the computational cost of the multi-frontal direct solver algorithm executed over such tensor product grids. We present an algorithm for estimation of the number of floating-point operations per mesh node resulting from the execution of the multi-frontal solver algorithm with the ordering obtained from the element partition trees. Next, we propose an algorithm that introduces C0 separators between patches of elements of a given size based on the stimated number of flops per node. We show that the computational cost of the multi-frontal solver algorithm executed over the computational grids with C0 separators introduced is around one or two orders of magnitude lower, while the approximability of the functional space is improved. We show O(NlogN) computational complexity of the heuristic algorithm proposing the introduction of the C0 separators between the patches of elements, reducing the computational cost of the multi-frontal solver algorithm

Petri Nets Modeling of Dead-End Refinement Problems in a 3D Anisotropic hp-Adaptive Finite Element Method

We consider two graph grammar based Petri nets models for anisotropic refinements of three dimensional hexahedral grids. The first one detects possible dead-end problems during the graph grammar based anisotropic refinements of the mesh. The second one employs an enhanced graph grammar model that is actually dead-end free. We apply the resulting algorithm to the simulation of resistivity logging measurements for estimating the location of underground oil and/or gas formations. The graph grammar based Petri net models allow to fix the self-adaptive mesh refinement algorithm and finish the adaptive computations with the required accuracy needed by the numerical solution